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ABSTRACT 

We present the results of an intensive photometric and spectroscopic monitoring campaign 
of the WN4 Wolf-Rayet (WR) star WR1=HD 4004. Our broadband V photometry covering a 
timespan of 91 days shows variability with a period of P=16.9^Q 3 days. The same period is also 
found in our spectral data. 

The light-curve is non-sinusoidal with hints of a gradual change in its shape as a function of 
time. The photometric variations nevertheless remain coherent over several cycles and we esti- 
mate that the coherence timescale of the light-curve is of the order of 60 days. The spectroscopy 
shows large-scale line-profile variability which can be interpreted as excess emission peaks moving 
from one side of the profile to the other on a timescale of several days. 

Although we cannot unequivocally exclude the unlikely possibility that WR 1 is a binary, we 
propose that the nature of the variability we have found strongly suggests that it is due to the 
presence in the wind of the WR star of large-scale structures, most likely Co-rotating Interaction 
Regions (CIRs), which are predicted to arise in inherently unstable radiatively driven winds when 
they are perturbed at their base. We also suggest that variability observed in WR6, WR134 
and WR 137 is of the same nature. Finally, assuming that the period of CIRs is related to the 
rotational period, we estimate the rotation rate of the four stars for which sufficient monitoring 
has been carried out; i.e. v rot =6.5, 40, 70 and 275 km s" 1 for WR1, WR6, WR134 and WR137, 
respectively. 

Subject headings: stars: Wolf-Rayet — stars: individual(WR 1, HD 4004) — stars: winds, outflows — 
stars: rotation 
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1. Introduction 

The winds of Wolf-Rayet (WR) stars are 
well-known to be non-uniform on small phys- 
ical scales. Because of the inherently unsta- 
ble nature of radiatively-driven ho t stellar winds 



(|Owocki. Castor fc Rvbicki 1988). instabilities, 



which reveal themselves in the spectrum as narrow 
excess emission peaks superposed on the broad 
emission lines, appear stochastically, propagate in 
the wind and disap pear after several hours (e.g. 



Moffat et al. 1 119881 ). This dumpiness leads to 



emission lines that ar e variable at a level of up to 
- 5 % of the line flux (jSt-Louis et al.ll2009h . These 
changes are, of course, random and no periodicity 
is expected. 

In certain cases, WR stars also display large- 
scale line-profile variability (lpv). For some stars 
this spectroscopic variability has been shown to be 
periodic and to originate from a different physical 
processes. Of course, known massive WR+O bi- 
naries produce clear periodic radial-velocity (RV) 
variations from the orbital motion of the stars. It 
has also been realized in the past two decades that 
in such systems the winds from the stars collide, 
forming a shock cone that wraps around the star 
with the smaller momentum flux (generally the 
O star). This interaction induces kinematically- 
characteristic variations in the line profiles that 
are also periodic. 

However, large-scale lpv has also been found 
in WR stars that are not known to be binaries. 
The two most studied cases, WR6 (WN4) and 
WR134 (WN6), have been monitor ed intensively 



i n photometry and spectroscopy by iMorel et al. 



(|l997t Il998l Il999ah . In both cases, these authors 
have observed periodic large-scale lpv and a com- 
plex light-curve with the same periodicity. To 
that list of two stars, we may also add WR137, 
a WC7pd star in a long-period binary system. 
This star shows large-scale lpv with a period of 
1.2 d that is unlikely to be related to its 09 
companion nor to the wind-wind collision zone 
( Lefevre et al. I l2005al ) . This behaviour is most 
likely explained by the presence of large-scale den- 
sity structures in the wind, such as co-rotating 
inter a ction regions (CIRs) ( Cranmer k, Owockil 



1996; iFullerton et al. 1 119971) . but the possibility 



of the presence of a compact companion or of 
a low-mass main-sequence star cannot be com- 



pletely excluded. 

If the large-scale lpv and the photometric vari- 
ability of single WR stars can be associated with 
CIRs, it may imply that the period of these vari- 
abilities corresponds directl y to the rotational 
period of the star. Indeed, ICranmer fc Owocki 
(|l996j ) propose that "spots" fixed to the stellar 
surface, caused either by pulsations or magnetic 
field activity, are at the origin of the CIRs. Hence, 
once the radius at which the CIR originates is 
known, we can determine the rotational velocity 
of the star at that point. It follows that the ro- 
tation rate of at least some WR stars could in 
principle be determined by carrying out a system- 
atic investigation of the v ariability of all single 



St-Louis et ail (|2009l) 



WR stars. As a first step, 
and Chene & St-Louis (in prep) set out to identify 
new candidates for CIR-type variability by con- 
ducting a survey of all apparently single Galactic 
WR stars brighter than v ^12. 5 th magnitude. For 
each star in their sample, they obtained 4-5 spec- 
tra which allowed them to establish a list of WR 
stars showing large-scale lpv. The next step is to 
observe intensively each of the candidates in order 
to verify and determine the periodicity. 

In the above-mentioned survey, WR 1 was one 
of the most striking cases of large-scale lpv with 
changes reaching 8-10% of the line flux and easily 
distinguishable large scale subpeaks superposed on 
the broad wind emission profiles. Consequently, it 
was the first WR star on which we concentrated 
our efforts. Several previous studies claiming un- 
reproducible periods for th e variability of this sta r 
have been summarized 
Later 



Nicdzielski 



Morel et al. I (|1999bh 



2000al lbl) demonstrated that 
the spectrum of WR 1 varied greatly from night- 
to-night reaching 50% in the equivalent width of 
the HeiA5876 line. On the other hand, they found 
that line-profile changes during one night are much 
smaller with a typical scatter of 0.5 A in equivalent 
width for the HeiiA5411, CivA5808 and HeiA5876 
lines. The changes, however, were found to be sys- 
tematic during the course of an entire night and 
reach a total of 3-4 A. Moreover, the search for pe- 
riods smaller than 2 days failed and only indica- 
tions of long-term variability could be suggested. 
The investigation of photometric va riability over 



more t han 16 days carried out by IMorel et al. 



(|l999bl ) did not lead to the identification of a pe- 
riod. Consequently, the period search must be 
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done using data taken over a time range of at least 



Q2007D 



twice as long. More recently, iFlores et al 
claimed a period of P=7.684 days in lpv of this 
star which they attribute to the ejecta of streams 
or jets from the stellar surface. 

In this paper, we present the results of an in- 
tensive monitoring campaign of WR 1 extending 
over several weeks in photometry and spectroscopy 
with the aim of determining the nature of the 
variability and, eventually, if found to be associ- 
ated with CIRs as suspected from our survey ob- 
servations, the rotation period of WR1. In Sec- 
tion [2j we present our photometric and spectro- 
scopic observations and the data reduction proce- 
dures. In Sections [3] and H] we describe our results, 
and in Section [5] we discuss the possible interpre- 
tations. Finally, our conclusions are presented in 
Section 15.31 



2. Observations 
2.1. Photometry 

We monitored WR 1 in broadband V using 
CCD-imagery at the 0.81m Tenagra Observatory 
ltd. from 20 November 2003 to 18 February 2004. 
During this period covering 91 nights, three frames 
were obtained in succession every clear night with 
an exposure time of 60s. Three additional frames 
with a 45s integration time were obtained at an 
airmass of ^1.45 during the first 43 nights and 
three others with a 30s integration time at an air- 
mass of ~1.65 during the first 11 nights. Due to 
problems with the detector in December 2003- 
January 2004, which resulted in a gap in our 
dataset of almost 20 days, 33 frames were unus- 
able. Also, 87 others were rejected due to bad 
seeing and/or low transparency of the sky, giving 
finally a total of 231 usable frames. In order to 
correct for the differential refraction of the atmo- 
sphere, we observed three frames in broadband B 
for each of two very clear nights at Tenagra Obser- 
vatory ltd., i.e. the 6 th and the 17*' 1 of December 
2003, for a total of 6 frames. 

All images were uniformly reduced using stan- 
dard procedures carried out with routines writ- 
ten in Interactive Data Language (IDL). After the 
bias, dark and flat-field treatments, we performed 
aperture photometry using the APER Astrolib rou- 
tine in IDL on all stars present in the field of WR 1. 
We adopted an aperture size equal to twice the 



FWHM of the point spread function (PSF) and an 
annulus of sky was selected with inner and outer 
radii of respectively 4 and 8 times the FWHM of 
the PSF. 

Once the photometry was obtained for all 
frames, we then applied a correction for the differ- 
ential refraction of the atmosphere. To do so, we 
needed to know how the flux of a star varies with 
airmass and according to the slope of its spec- 
trum. For each star in our field of view, a relative 
(B — V) color was obtained by subtracting the av- 
erage fluxes measured in the B and U-band frames 
observed at the same airmass. We then obtained 
the slope M of the linear relation describing the 
variability of the flux as a function of airmass for 
each star during the best clear nights. Although a 
certain number of stars can present intrinsic flux 
variations within a single night, we can assume 
that it is not the case for most of them. Using a 
robust least absolute deviation method, we fitted 
M( B _y), the slope of the flux-airmass relation ap- 
propriate for a given (B-V) color, and obtained a 
linear expression for the correction that must be 
applied to the flux as a function of airmass. The 
final corrections are contained in a range of 0.1 to 
1 %. 

The final light-curves were obtained by sub- 
tracting the magnitude of WR1 to that of two 
comparison stars (see Figured]). Since no star with 
a magnitude similar to that of WR 1 was avail- 
able in the field of view, we decided to average the 
brightest stars showing small light deviations. Ten 
stars were selected and divided in two groups to 
create two comparison stars of similar magnitudes. 
The accuracy of the final light-curve of WR 1 var- 
ied during our observing period. We have thus 
estimated a typical error for each day as follows : 
for a time interval of 10 days centered on a given 
day of observation, we have calculated the stan- 
dard deviation (a) of the cl-c2 light curve. All 
values are shown in the bottom panel of Figure [1] 
For the entire observing period, the standard de- 
viation is cr=0.012 mag. 
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Tabic 1: Spectroscopic observing campaigns 



Run # 


Telescope 


Dates 


A Coverage 


AA 


No. of sp. 


SNR 






(UT) 


A 


(3 pix.) 






1 


OHP 1.52m 


2003/07/07-14 


5230-6140 


1.4 A 


30 


135 


2 


CFHT 


2003/08/08-13 


4380-6240 


1.5 A 


33 


185 


3 


OMM 


2003/08/22-24 


4900-6020 


1.6 A 


9 


130 


4 


OMM 


2004/08/21 - 


4310-6440 


1.6 A 


170 


145 






2004/09/06 




1.4 A 






5 


DAO 1.85m 


2004/09/18 - 


5140-5980 


84 


90 






2004/10/13 
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2.2. Spectroscopy 

We collected a total of 326 spectra of WR1 
during five dedicated runs distributed in July- 
August 2003 and August-October 2004 with the 
1.52m telescope of the Obervatoire de Haute- 
Provence (OHP), the 3.6m Canada- France- Hawaii 
Telescope (CFHT), the 1.6m telescope of the Ob- 
servatoire du Mont Megantic (OMM) and the 
1.85m telescope of the Dominion Astronomical 
Observatory (DAO). The details of these observ- 
ing runs are summarized in Table [1] We list the 
run number, the telescope used, the dates of the 
observations, the wavelength coverage, the spec- 
tral resolution, the number of spectra obtained 
and the average signal-to- noise ratio (SNR). Un- 
fortunately, no spectra were obtained simultane- 
ously with our photometry described in the previ- 
ous section. 

The bias subtraction, flat-fielding, spectrum ex- 
traction, sky subtraction and wavelength calibra- 
tion of all spectra were executed in the usual way 
using the irai£] software. Calibration-lamp spec- 
tra were taken every 30-40 mins depending on 
the run. The accuracy of the wavelength calibra- 
tion estimated by measuring the wavelength of 10 
lamp emission lines is variable depending on the 
instrument used. It is 0.005 A for OHP, 0.2 A for 
A < 5200 A and 0.02 A for A > 5200 A for CFHT, 

0. 01 A for OMM and 0.04 A for DAO. Special 
care was taken for the normalization of the spec- 
tra. First, a mean was made for each run. Then 
each spectrum of a run was divideded by the run 
mean and the ratio fitted with a low order Legen- 
dre polynomial (between 4 th and 8 th order). The 
original individual spectrum was divided by this 
fit, and was therefore at the same level as the run 
mean. When this procedure was done for each 
run, the run means were then put at the same 
level by using the same procedure as described 
above. This allowed us to put all individual spec- 
tra at the same level. Then, we combined all run 
means into a global, high quality mean, which was 
then fitted in selected pseudo-continuum regions, 

1. e. wavelength regions where large emission-lines 
do not dominate. More specifically, these regions 



are : 4376.0 - 4381.0A, 4473.0 - 4477.0 A, 4798.0 

- 4816.0 A, 5002.0 - 5147.0 A, 5349.0 - 5362.0 A, 
5558.0 - 5663.0 A, 5956.0 - 5992.0 A and 6295.0 

- 6321.0 A. These regions are shown in Figure 21 
which we will discuss later in Section 2J Finally, 
the fitted continuum function is applied to each 
individual spectrum. The error on the continuum 
normalization measured as the standard deviation 
of individual spectra around the continuum func- 
tion is typically of 0.5%. 

3. Photometric Variations 

The light-curve of WR 1 is plotted as a func- 
tion of the Heliocentric Julian Date (HJD) in 
Figure [1] The flux from the star alternatively 
increases and decreases with an amplitude that 
varies between AV=0.06-0.12 mag within 5-7 days 
(<r(cl — c2) ~ 0.012 mag), suggesting that the 
changes are periodic with a period of at least 
5 days. The variability within a single night 
is small and typically reaches an amplitude not 
higher than 2a. Since only slow variability is 
observed in pho tometry and spectro s copy within 



a sing le night (jMorel et al. I Il999bt iNiedzielski I 
2000af ). we binned all the data obtained during 
a given night to increase the SNR. 

We performed a period search using two inde- 
pendent methods. The first is a ph ase-dispersion 
minimization (PDM) algorithm ( Stellingwerf I 



19781 ) which is well suited to cases in which only 
a small number of observations are available over 
a limited period of time, especially if the signal 
is highly non-sinusoidal. The method consists 
in folding the data points in phase using differ- 
ent trial periods and dividing the resulting curves 
in a predetermined number of phase bins. For 
each trial period we can define Sj, the variance 
of bin j, S 2 , the variance of all S 2 and finally 
O = S 2 /(J 2 , where a 2 is the variance of all data 
points. The other method used is the CLEA N algo- 



rithm (jScargle 1119821 : iRoberts et al. 1119871 ) which 



1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories (NOAO), which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) 
under cooperative agreement with the National Science 
Foundation (NSF). 



has the advantage of taking into account the un- 
evenness of the data sampling, since it "cleans" 
the discrete Fourier spectrum with a window func- 
tion. Both, the O spectra from the PDM method 
and the power spectra from the clean method 
are shown in Figure [2] for WR 1-cl, WR l-c2 and 
cl-c2. 

The total time elapsed between the first and 
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Fig. 1. — Light-curve of WR1. The top and mid- 
dle panels show the difference between the mag- 
nitude of WR 1 and that of the "artificial" com- 
parison stars cl and c2 respectively as a function 
of time. Both cf and c2 were constructed from 
the average of five bright stars. The 6 cycles of 
the 16.9-day period covered by the photometric 
dataset are delimited. The bottom panel shows the 
difference between the magnitude of cl and c2 as 
a function of time. A series of vertical lines show 
the typical error as a function of time (see text for 
details). 

last point of the light-curve is 91 days and, as men- 
tioned above, we combined the data to obtain one 
measurement per night. Therefore, the search was 
performed over a range of frequencies from 2/91 to 
0.5 day -1 with a step of 0.00111 day -1 , respect- 
ing the Nyquist criterion. When using the PDM 
method for (WR 1-cl) and (WR l-c2), one fre- 
quency (vi,pdm) is found to be predominant over 
all others. Interestingly, when only the first 60 
nights are included in the period search, vx,pdm 
becomes even more significant, and 3 other fre- 
quencies can be distinguished from the noise in 
the periodogram. The strongest detection is at 
^i,P25M=0.059io!oo2 day -1 and the three others 
are at V2,PDM=0.028±° f Q l is 3 ,pdm =0.118±g;ggJ 
and ^4,pum=0.236±o:oo2 day -1 . Note that the 
peaks at V2,pdm, vs,pdm and va,pdm reach com- 
parable levels and are much smaller than the 
peak at v\,pdm- The periodogram of (cl-c2) 
shows no signal at all, indicating that there is 
no periodic variability present in the light-curve 
of either of the comparison stars. The power 
spectra obtained using the clean algorithm for 



Fig. 2. — Periodograms: a) 9-spectrum obtained 
from the PDM analysis of the WR1 - cl curve. 
The four frequencies Vi,pdm, vi,pdm, vs,pdm 
and va.pdm are indicated with an arrow, b) Same 
as a) for WR l-c2. c) Same as a) for cl-c2. d) 
CLEANed power spectrum of WR 1-cl (green), 
WR l-c2 (bleu) and cl-c2 (red). Three levels of 
confidence at 99%, 95% and 90% are plotted in 
black. The three frequencies v\,clean, ^2,clean 
and vzfiLEAN are indicated with an arrow. 

both the (WR 1-cl) and (WR l-c2) light-curves 
show only three peaks that reach values higher 
than the 99% confidence level threshold and that 
are absent in the comparison light-curve (cl- 
c2). These peaks are at z/ 1) c i z,.E!>iAr=0.119±0.007 
(~ V3,pdm), ^ 2: cl£Mat=0.057±0.007 (~ vi, PD m) 
and z/ 3iC . L£Mi v=0.240±0.007 day -1 (~ v^ PDM ). 

The two main frequencies identified using the 
two different methods (^=0.059, 0.119 day -1 ) can 
be linked by an integer number. It is therefore 
clear that one is an harmonic of the other. Ac- 
cording to the PDM analysis, the frequency with 
the strongest signal is Vi,pdm=0-Q59_q]oo2 day -1 
(period of 16. 9^3 days), but according to the 
CLEANed spectrum, it is ^i i clsatv=0.119±0.007 
day -1 (period of 8.4±0.5 days). In an attempt to 
settle which among these two possible periods for 
WR 1 is the correct one, we folded the light-curve 
in phase using the two different frequencies. The 
frequency that produced the smallest dispersion 
and the clearest curve is vx,pdm- This result is 
presented in Figure [3] where each 16.9-day cycle of 
(WR l-c2) is plotted using a different symbol. The 
6 cycles of 16.9 days covered by our photometric 
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dataset are delimited in Figure [TJ The 16.9-day 
folded light-curve for cycles 1 to 4 shown in the 
top panel of Figure [3] shows basically three bumps; 
one broad and small centered at cf> ~0.6 and two 
stronger ones spaced closer in phase at <j> ^0.95 
and 0.15. The bottom panel shows the same data 
folded with a 8.4-day period. Clearly the scatter 
is much larger than when folded with the shorter 
period. A careful examination of Figure [T] also 
shows that there is no clear pattern repeating in 
any successive 8.4 days, which supports this re- 
sult. The fifth and sixth cycles demonstrates the 
epoch-dependancy of the light curve. Indeed, the 
curve shown in the middle panel of Figure [3] shows 
the disappearance of the second bump near phase 
0.95. Indeed, the third bump seems significantly 
higher and wider than during the previous cy- 
cles and starts from a lower level, but the star 
seems to have skipped the previous increase and 
decrease in flux, as if the two nearby peaks had 
merged. We have no data at phase 0.6 for these 
cycles so we cannot tell if the bump at that lo- 
cation is still present. To show the difference be- 
tween the shapes of the light curves, we overplot 
in the middle panel showing cycles 5 and 6, the 
cycle 1 to 4 light curve as a thick grey curve which 
roughly emcompasses all data points. This epoch 
dependency of the light curve of WR 1 shows that 
a period can only be found when observations are 
taken consecutively, during a period of time for 
which the physical conditions of the region where 
the continuum flux originates do not change. This 
means that no data taken after a certain coherence 
timescale can be added to the sequence in order 
for a period search to be successful. Indeed, even 
if the period remains the same, the amplitude and 
the shape of the light-curve will change. We can- 
not determine with accuracy the time of coherence 
for WR1, but according to our data, it seems to 
be at least 4 cycles, i.e. ^60 days. 

A li ght- curve was obtained by iMorel et~al 



(|1999bT ). During their observations extending over 
a period of 16 days, the authors found only one 
bump lasting 5-days in width, followed by a pe- 
riod of 11 days during which the curve was flat. 
In the context of an epoch-dendendant light curve 
shape, it is interesting to note that their curve is 
not inconsistent with our choice of adopting a pe- 
rio d of P~ 16.9 days as the single bump observed 



width to the bump obtained in the present data 
during cycle 6 near phase 0.15. This is shown in 
the middle panel of Figure |3] where we have arbi- 
trary shifted the curve in phase. On the other 
hand, the fact that the star remained constant 
during 11 days does not bring support to a shorter 
8.4-day period. 



• cycle 1 * cycle 3 
' cycle 2 ■ cycle 4 



> cycle 6 + Mor99 



phase (16.9 days) 



_ a -k o 



by IMorel et al. (1999b) is similar in height and 



phase (8.4 days) 

Fig. 3. — Top panel: V-band light-curve of WR 1 
folded with the 16.9 days period. Phase is ar- 
bitrary. The first cycle is represented by circles, 
the second by stars, the third by triangles and the 
fourth by squares. Middle panel: The light-curve 
from the top panel is reproduced here and drawn 
with a continuous thick, gray line. Overplotted 
to it are the data from cycles 5 (triangles) and 
6 (diamo nds). The na r rowban d v light-curve ob- 
tained bv lMorel et al. ( 1999bl ) is also plotted with 
+ signs. Bottom panel: The same data folded with 
a 8.4-day period. The symbols represent the 9 dif- 
ferent cycles of the 8.4-day period. 

If the period of 16.9 days is the real one, why 
did the CLEAN analysis yield a period with half this 
value? The answer may reside in the way the algo- 
rithm works. Indeed, the discrete Fourier analysis 
looks for a sinusoidal signal. Now, the light-curve 
of WR 1 we observed can be interpreted as two 
nearly similar sinusoidal curves per cycle. This is 
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the reason why the CLEANed spectrum has more 
power in the peak at v\,pdm /2. However, when 
this frequency is used to fold the light-curve in 
phase, the smaller bump is superposed on the big- 
ger one, and a large dispersion around the curve 
is found. This larger dispersion explains why the 
PDM analysis considered the V\,pdm/1 result of 
lower interest. From this and mor e importantly 
when t aking into account the data of lMorel et al. 



the S spectrum as follows : 



(|1999b[ ) we conclude that the most probable inter- 



pretation is that the period of 16. days, ob- 
tained using both methods, is the real period and 
that the other detections are simple harmonics or 
a consequence of the fact that the light-curve is not 
a simple si nusoide. Note that the 7.684-d period 
claimed bv lFlores et al. ( 2007 ) has been obtained 
from spectra taken during two different runs of 5 
and 8 nights separated by ~6 years. Due to an 
insufficient time coverage and to the complexity 
of the variability pattern, they have detected only 
a fraction of the real period. 

4. Spectroscopy 

The photometry discussed in the previous sec- 
tion provides information on the variability of the 
continuum of the star since for this star only 10% 
of the flux in broadband V comes from emission 
lines. This continuum originates in deep regions 
of the stellar wind. In order to probe the wind 
at different radii, we must study spectral variabil- 
ity, because the different WR emission-lines are 
formed in regions of different velocity regimes (i.e. 
different radii), depending on the io nisation po- 



tential of th e given elements (e.g. iKuhi 1973: 



Hillierlll987h . From spectral variability, we hope 



to obtain clues as to the origin of the period found 
in photometry. 

4.1. Level of Variability 

In Figure 2^,, we present the mean spectrum 
for each of our spectroscopic data runs listed in 
Table [T] (thick black lines). Line identification is 
provided in the panel for run 1. 

To illustrate which parts of our spectra are vari- 
able, we have calculated the temporal variance 
spectru m (TVS) according to t he method intro- 
duced by Fullerton et al. I ( 19961 ). Then, assuming 
that our data are governed by a reduced \ 2 distri- 
bution with N-l degrees of freedon, we calculated 



£,■(99%) = 



(TVS), 



2 2 



i(99%) 



(1) 



where N is the total number of spectra and do a 
standardized dispersion. The £ spectrum for each 
run is overplotted (in thin blue line) on the mean 
spectra in Figure A £(99%) value of 3, for 
example, means that we are 99% confident that 
that part of the spectrum is variable at a 3<r level 
(see ISt-Louis et al.l l2009l for more details). It 



can be easily seen that all lines show variability 
at high sigma levels with a high degree of confi- 
dence for each run. Note that since the TVS com- 
pares of the variability level at each pixel with the 
variability calculated in the continuum regions, no 
changes can be identified in the continuum with 
this method. 

In Figure B b, we present th e a sp ectrum which 
we defined in ISt-Louis et ail (|2009l ). This spec- 
trum gives the fraction of the line flux that is vari- 
able (and therefore is only defined within spectral 
lines). We present 4 regions of our wavelength 
range which include spectral lines. Note that each 
region is not covered by the same number of spec- 
tra. Hence, some are more noisy such as HeiiA4686 
for which we only have ~20 spectra as opposed to 
HeiiA5411 for which we have ^320. The reason 
we have so little spectra for the HeiiA4686 line is 
that we decided to allow this line to saturate in 
run 4, the run with the largest number of spectra, 
in order to increase the total number of spectra 
we obtained. Indeed, the alternative would have 
been to divide the exposure time by a factor 3 and 
observe three times as often, which would have led 
to tripling the amount of time spent reading the 
detector. It can be seen that all lines vary typ- 
ically at a ^7 % level. This is exac tly the same 
result found bv lSt-Louis et al. I (l2009h for this star, 
but based only on 5 spectra. 

4.2. Variability Pattern of Different Spec- 
tral Lines 

In order to compare the variability pattern of 
different emission lines, we have calculated cor- 
relation coefficients between each velocity bin of 
the different emission profiles, using the Spear- 
man rank-order correlation. This procedure yields 
a matrix of correlation that assesses how well the 



S 



variability pattern of two lines correlate. For lines 
with overlapping formation regions (which if often 
the case) , a correlation should be found if an over- 
density is present in the wind. In the case of line- 
formation zones far-removed from the acceleration 
region, large sections of the wind are associated 
with a narrow velocity interval. In such 
correlation would support "solid" rotating CIRs. 
In the case of perfectly correlated variations, a ma- 
trix unity is obtained. In Figure [SJ we present the 
correlation matrices of HeiiA5411 with HenA4860, 
NvA4945 and CivA5808. The lowest level repre- 
sents a significant correlation at the 99.5% confi- 
dence level. The correlation matrix of HeiiA5411 
with HeiiA4686 is not shown he r e, but has already 
been published bv lMorel et al. I (|l999bl ). The vari- 
ability pattern of HciiA5411 is clearly correlated 
with that of all other emission lines, except for 
HeiA5876 and for NvA4603/20, which is blended 
with the highly variable absorption part of the 
P Cygni pro file of HeiiA46 8 6. One can note that, 
contrarily to iMorel et al. 1 (|l999bh . we were able 
to detect a clear correlation between the variabil- 
ity pattern of the HeiiA5411 and NvA4945 lines. 
This is most likely due to the higher signal-to- 
noise ratio and the greater number of spectra in 
our dataset. 

In view of the above results, we only show here 
the variability pattern of the HenA5411 line which 
is a strong and isolated line and is common to all 
our observing runs. The variability of that line 
when observed during individual nig hts is typical 
of what has already been reported by iNiedzielski I 
(2000a) for WR 1, i.e. low level variations occuring 
slowly and systematically during the night. How- 
ever, the changes are more interesting when dis- 
played in the form of a grayscale plot covering a 
time period of more than a week, as shown in Fig- 
ure El where we plot the complete dataset of our 
runs 4 and 5. For clarity, the spectra obtained dur- 
ing each night have been repeated once to fill the 
space in the grayscale plot corresponding to day- 
light, when no observations were possible. This 
way, no free space has been left between two con- 
secutive nights,which makes it easier to see the 
structures moving on a timescale of several days. 
Although our time coeverage is not ideal, one can 
begin to distinguish the typical "S-type" pattern 
caracteristic of CIR-type variability. 












A 







5800 5901 



Fig. 4. — a) Mean spectrum for each of our spec- 
troscopic data runs (thick black lines). The £ 
sp ectrum, calculate d with the method introduced 
by iFullerton et al~l (|l996l ) and described in text, 
is overplotted for each run (thin blue lines). Line 
identification is provided in the panel for run 1 
and the regions adopted for the rectification of 
the spectra are shown in t he panel for ru n 5. b) 



(2009) 



a spectrum as defined in ISt-Louis et al. 
calculated in the wavelength regions where the 
strongest emission lines are present. The dotted 
line represent the mean a-value for all emission 
lines across the observed spectrum. 

4.3. Radial Velocity, EW, Skewness and 
Kurtosis Variations 

We next investigate the variability of the RV of 
spectral lines as well as the changes in the global 
intensity of the line profile by measuring the equiv- 
alent width (EW), the skewness and the kurtosis 
of the HeiiA5411 line. The radial velocities were 
obtained by cross correlating individual spectra in 
the wavelength interval AA=5140-5980 A with the 
mean spectra. The RVs of the different campaigns 
were obtained separately. We corrected the sys- 
tematic shifts between runs by cross-correlating 
the mean spectra of the runs. The typical val- 
ues for the applied corrections ranged between - 



9 



HellA4860 NVA4945 CIVA5808 




-1500 1500-1500 1500-1500 1500 
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Fig. 5. — Contour maps of the correlation ma- 
trices of HeiiA5411 with HeiiA4860, NvA4945 and 
CIVA5808. The line contours indicate a correla- 
tion (in black) or an anti-correlation (in red) in 
the variability pattern of the lines. The lowest 
contour level represents a correlation at the 99.5% 
confidence level. 

20 km s -1 and -40 km s -1 . The EWs of emis- 
sion lines were calculated by integrating the func- 
tion (1 — F\), where Fx is the rectified line flux 
between AA=5348-5460 A. In order to highlight 
the variability, we prefer to use AEWs obtained 
by dividing the EWs by the mean value of all 
EWs. In order to verify that the AEW of the 
HeiiA5411 line is representative of the curve for 
any other emission line in the wavelength range 
of our observations, we compared its AEW to the 
one of the NvAA4603/20, HeiiA4686, HeiiA4860, 
CivA5808 and HeiA5876 lines. We present in Fig- 
ure [7] the AEW of the HeiiA4686 and HeiiA5411 
lines observed in run 2. The changes clearly fol- 
low the same pattern. In general, except for 
the blended NvAA4603/20 and the weakly varing 
HeiA5876 lines, with the quality of the data in 
hand, there is no significant difference that can be 
seen in the AEW value of all observed lines. The 
other moments of the HeiiA5411 line are calculated 
in the wavelength interval AA=5348-5460 A. The 
n th central moment is defined as follows: 

H n = Z j (\ j -\) n I j /X j I j (2) 

where 

A = ZjXjlj/Zjlj, (3) 

with Ij being the intensity of the line and Aj the 
wavelength. The skewness is ^3/ ' p^ 2 and the kur- 
tosis is ^4/^2- All these quantities are plotted as 
a function of time in Figure [51 

On average, the EW of the HeiiA5411 line in- 
creased by ~15 % during the first three runs (left 




-3000 3000 

velocity (km/s) 



Fig. 6. — Top panel: Grayscale of residuals for the 
HeiiA5411 line obtained by subtracting the mean 
spectrum of all runs from the rectified spectra 
taken during runs 4 and 5. Bottom panel: Super- 
position of all spectra taken during the two runs. 
Note that for clarity, the time corresponding to 
daylight has been filled by the duplicate of the 
spectra obtained during the previous night. 
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Fig. 7.— AEW of the HeiiA4686 (diamonds) and 
HeiiA5411 (triangles) lines observed in run 2. The 
error bars are 3 a long. Note that due to probems 
with the first spectrum of that run at A <5200A, 
AEW could not be computed for HeiiA4686 on the 
first night. 

panel) and by ^5 % during the last two (right 
panel). The RVs show long-term changes with a 
full amplitude K^70 km s" 1 as well as shorter 
timescale changes which are particularly appar- 
ent during run 4. However, the skewness, which 
characterizes the asymmetry of the line profile, 
is strongly anti-correlated with the RVs. Indeed, 
these values are anti-correlated with a 90 % confi- 
dence level when we remove run 5, which seems to 
have abnormally high dispersion. This strong anti- 
correlation indicates that the RV changes that we 
measure essentially come from changes in the de- 
gree of asymmetry of the line profile of HeiiA5411, 
which dominates the spectrum in the wavelength 
interval in which the cross-correlation was carried 
out. Consequently, if WR 1 shows intrinsic RV 
variability in HeiiA5411, it must be at a very small 
level and therefore is difficult to detect in view of 
the line-profile variability. Both EW and RV of the 
HeiiA5411 line vary with an amplitude of at most 
~5 % in a single night. This can b e comp ared 
with the observations of iNiedzielski (2000a) who 
found EW variations of HeiiA5411 of ~15% for 
their first run with a structured variability pat- 
tern over a At=5 day period, but only for 
their second run on a timescale of ~2 days. Un- 



fortunately, iFlores et al. I (|2007l ) did not observed 
this line. The variations in EW are clearly epoch- 
dependant. 
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Fig. 8.— Top panel: AEW of the HeiiA5411 
line, obtained by integrating the line flux and di- 
viding the result by the mean value of all the 
EWs. Second panel from top: RVs obtained by 
cross-correlating the spectra in the 5140-5980 A 
range. Second panel from bottom: Skewness of the 
HeiiA5411 line calculated as skewness=^3/'/Z2 / ' 2 , 
where fj, n is the central moment of order n of the 
line. Bottom panel: kurtosis of the HeiiA5411 line 
calculated as kurtosis=/Z4//z|. 



4.4. Search for Periods 

A period search was performed using the PDM 
method on the EWs, skewness and kurtosis vari- 
ations (but not on the RVs, due to these simi- 
larity with the skewness). No clear period could 
be found in the EW curves. However, the peri- 
odograms calculated for both the skewness and 
the kurtosis curves have a significant peak at 
0.059 day~\ i.e. the photometric frequency 
vi.pdm discussed in the previous section (see Fig- 
ure |9]). The curves folded with that frequency 
are shown in Figure 1101 Each point is the mean 
value for a single night. Run 1 is represented by 
circles, run 2 by squares, run 3 by triangles, run 

4 by diamonds and run 5 by stars. The folded 
curves for the AEW and the skewness have quite 
a large scatter and poorly confirm the detected 
period. However, the folded curve for the kurtosis 
is more promising. The largest scatter is found 
between phases 0.55 and 1.0, but it can be greatly 
reduced if one conciders separately runs 1, 2 and 
3 observed in 2003 (At ~ 50 days) and runs 4 and 

5 observed in 2004 (~ 55 days). Interestingly, this 
is also the time interval aftehr which a change in 
the shape of our lightcurve was observed. One 
maximum, centered at phase^O.l is common to 
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all runs. The variations between phases 0.55 and 
1.0 can be described as one additionnal maximum 
for runs 1, 2 and 3 and a minimum (around phase 
0.6) followed by a maximum (around phase 0.9) 
for runs 4 and 5. 




frequency (days ') 



Fig. 9. — 8-spectrum obtained from the PDM 
analysis of the skewness curve (black, dashed 
line) and the kurtosis curve (red, solid line) of 
HeiiA5411, presented in Figure [H The photomet- 
ric frequency v\.pdm discussed in Section [3] is in- 
dicated with an arrow. 

It is not necessarily trivial to link the variability 
of the global moments of the lines and the wind 
variability of the star. Here, we present an in- 
vestigation that aims to reveal the periodicity of 
the movement of the extra bumps observed super- 
posed on the emission line profiles. These can be 
seen moving from one side of the profile to the 
other for runs 4 and 5 in Figure [6] We define the 
scalar value a res , i.e. the standard deviation of 
the residuals obtained from the subtraction of two 
spectra observed separately with a time interval of 
At. Theoretically, when two spectra are obtained 
at the same point of a given variability pattern 
and therefore should be most similar, oy es reaches 
a minimum and takes the value of the quadratic 
sum of the noise levels of both spectra. Also, any 
significant difference between the two spectra will 
increase the value of a res proportionally to the 
magnitude of that difference. The great advan- 
tage of a res is that its value depends weakly on the 
epoch depedency of the signal. Indeed, even if the 
signal changes after 4 or 5 cycles (see Section [3]), 
all spectra from any cycle can be used at a given 
At, since only the spectra separated by the time 
interval At are used. The effect of the epoch de- 
pendency becomes significant for At greater than 
the coherence time, i.e. in this case, when two 
spectra separated by more than (at least) ~4 cy- 
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Fig. 10. — The AEW, skewness and kurtosis 
curves for HeiiA5411 folded with the frequency 
^i,pum=0.059 days -1 . Each point is the mean 
value for a single night. Run 1 is represented by 
circles, run 2 by squares, run 3 by triangles, run 4 
by diamonds and run 5 by stars. 

cles (^67 days) are compared. 

In Figure[TT]we plot the values of a res measured 
within a wavelength interval of 65 A centered on 
A=541lA for At < 55 days, using all the spec- 
tra from our 5 observing runs. The wavelength 
interval has been chosen to include the part of the 
HeiiA5411 line where most of the lpv occurs, with- 
out adding too much of the less intense parts that 
do not contain much information and have a much 
lower SNR. For clarity, we have organized the val- 
ues of a res as a function of At in a grayscale. 
In black are plotted the regions with the higher 
sumber of data points and, in white, the regions 
with no data point. Since the spectrum of WR 1 
does not vary very much within a single night (see 
above), we assume that the mean value of a res 
measured for all At < 1 day can be used to de- 
termine the minimum value for a res (plotted with 
a dotted-line) . Also, the spectra have been nor- 
malized in amplitude to avoid any change coming 
from the line dilution caused by the variation in 
the continuum that we observe in photometry. 

The periodogram calculated for the a res vs At 
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curve using the PDM method is shown in Fig- 
ure [12] The four frequencies v\,pdm, V2,pdm, 
v?,,pdm and v^pbm are indicated by arrows. In 
Figure [TTJ a vertical dashed-line is plotted at a 
constant interval of 16.9 days. Interestingly, a res 
reaches a minimum value near At=nxl6.9 days, 
where n=l, 2 and 3. These minima are signifi- 
cant when compared with a\ n , but they do not 
reach the minimum value corresponding to the 
noise level. This may indicate that even after one 
cycle, the variability pattern has changed slightly 
(but not drastically). Two other small minima 
can be seen at At ^4.5 and 11.8 days, but they do 
not repeat at greater At. This can be caused by 
line-profile variability that shows a similar pattern 
occuring more than once within the period for at 
least one cycle (as seen in the photometry). 

Finaly, in order to illustrate the line-profile vari- 
ability pattern, we display in Figure fT3l the resid- 
uals obtained by subtracting the mean spectrum 
of all runs from the mean spectra of each night 
of runs 4 and 5 for the NvA4945 and HeiiA5411 
lines (note that the NVA4945 line was not ob- 
served during run 5). The location on the y-axis 
of the residual spectra is determined by the phase 
(P=16.9 days) at which the spectrum was taken. 
Different colors are assigned to residuals from dif- 
ferent cycles. 

One can see from Figure [13] that the residuals 
at nearby phases are not always strictly identical 
(as expected from the above analysis). However, 
some features seem to survive for several cycles. 
The clearest case is that of the bump that ap- 
pears at v ^+250 km s _1 at phase 0=0.40 of the 
first cycle (HJD-HJD(2000)=1702.54), that moves 
toward negative velocities during phases 0=0.44- 
0.61 of the third cycle (HJD-HJD(2000)=1736.85- 
1739.83), and continues its blueward motion dur- 
ing phases 0=0.65-0.70 of the first cycle (HJD- 
HJD(2000)=1706. 70-1707.70) before disappear- 
ing. More moving structures can be followed on 
Figure [13] and are shown by dashed (for bumps) 
and dotted lines (for dips) . One can note that the 
motion of the bumps explains the changes we have 
observed in the central moments of the HeiiA5411 
line presented above. Indeed, the two biggest 
bumps trace a sinusoidal trajectory on the line, 
both in opposite directions. Hence, it is expected 
that the skewnesss, which measures the asymme- 
try of the line, would not vary too much due to 




Fig. 11. — The value of <j res measured for the 
HeiiA5411 line as a function of At. a res is defined 
as the standard deviation of the residual obtained 
from the subtraction of two spectra observed sep- 
arately with a time interval of At. The dotted line 
indicates the minimum value for a res due to the 
noise level. The dashed lines mark Ai=nxl6.9 
where n=l, 2 and 3. The vertical bar in the top 
left corner indicates the 1 sigma scatter of the a res 
value for a given At. 

the motion of the bumps, and that the kurtosis, 
which measures the degree of "peakedness" of the 
line, would be more sensitive to it. 

5. Discussion 



St-Louis et all (|2009h have already shown that 



WR l's spectrum shows large-scale lpv with a sim- 
ilar amplitude to what is observed for WR6 and 
WR 134, two presumably single WR stars showing 
periodic photometr i c and spectro scopic variability 



(|Morel et al.~lll997t Il998l Il999al lbh. With the in 



tensive monitoring campaign presented here, we 
are now able to conclude that, as for those two 
stars, WR 1 presents a unique period in its epoch- 
dependent light-curve and spectral varibility. The 
changes consist in large-scale bumps moving pe- 
riodically from one side of emission lines to the 
other. The origin of this type of variability is de- 
bated in several papers and the two main inter- 
pretations put forward are a binary system (with 
a compact or low-mass companion) and a rotat- 
ing non-spherically symmetric wind. In what fol- 
lows, we discuss these two possible interpretations 
in the context of the variability of WR 1 as char- 
acterized in this work, as well as that of WR 6 and 
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quency (days ') 



Fig. 12. — O-spectrum obtained from the PDM 
analysis of the o~ res curve. The photometric fre- 
quencies PDM) ^2, PDM 7 V3,PDM and V^ PDM 

discussed in Section [3] are indicated by arrows. 



WR 134 as detailed in the literature. We also dis- 
cuss the variability of the WR star member of the 
binary WR 137, which shows large-scale lpv with 
a 1.2d-period that cannot be caused by the 09 
companion nor the wind- wind collision zone, since 
the orbital period is ~13 years and the average 
distance between the two components is high. 

5.1. The Binarity Scenario 

Like WR6, WR134 and WR137, WR1 shows 
very small RV amplitude, consistent with null, 
and its soft X-ray luminosity is wi thin the nor- 



mal ra nge for single WR stars (see iMorel et al. 
1999b)) . Hence, it is completely excluded that its 



companion is an OB star. The presence of a com- 
pact companion, such as a neutron star or a black 
hole, in WR1, WR6 and WR134 was not see n 
as a viable in t erpret a tion bvllgnace et al. I (120031) , 



Morel et al. 




Skinner et al I (|2002al) ~and 



Morel et al. I (|1999af) . Indeed, the strong correla- 
tion between the lpv of Hen lines with that of a 
highly ionized line such as NvA4945 indicates that 
the ionizing shell around the compact compan- 
ion would have to be extremely large and would 
emit a high X-ray flux, which is not observed 



(IMorel et al. lll9 97. 1999a! IPollock" Habe rl fc Corcoran 
1995h . 

Therefore the only remaining possibility for 
the binary scenario for WR1, WR6 and WR134 
is that of a low-mass, non-compact companion. 
Indeed, assuming circular orbits and taking the 
published RV va riation amplitudes measu red for 



the three stars (jMorel et al. I Il997l Il998l and 



this work) as the maximum possible value (recall 
that RV variations are highly dominated by lpv) 







Fig. 13. — Residual of nightly mean spectra of 
runs 4 and 5 folded with a period of 16.9 days (left 
panel: NvA4945; right panel: HeilA5411). In black 
we plot the residuals from the first cycle covered, 
in blue the second, in green the third and in red 
the fourth. Moving structures can be followed and 
are shown by dashed (for bumps) and dotted lines 
(for dips). In the bottom panel, we plot individual 
intensity spectra superposed. 
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the periods lead to an upper limit for the mass 
of such a companion of 8 M© in all cases, if the 
inclination of the orbital plane is higher than 20°. 
As for WR 137, the scatter i n RV on a tim escale of 
a few days is ~10 km s _1 (jLefevre et al 



2005aD 



and its measured M sin 3 i=3.4±1.0 M©. Hence, 
if the inclination of the orbital plane is higher 
than 20° , the mass of the companion would have 
to be as small as 0.8 M in an orbit of 7.6 R©, 
which is very unlikely. But, can a low-mass com- 
panion be responsible for such a complex and 
epoch-dependant patte rn of photometric and spec 



trosco pic va riability? iMoreno fc K oenigsbcrger 



1999 ) and IMoreno. Koenigsberger fc Toledano 



2005) have shown that tidal interactions involv- 



ing a relatively low-mass companion in a binary 
system can produce very small-scale surface oscil- 
lations leading to lpv in photospheric absorption 
features. However, it is unclear how such a pro- 
cess can affect the massive wind of a WR star 
and lead to large-scale variability of the strong 
emission lines. It is possible that it could serve 
as a seed mechanism for CIRs. However, in such 
a case, it is likely that the period that would be 
detected would be the rotation period of the WR 
star rather than the orbital period of the binary, 
if not a combination of both. 

Therefore, although very unlikely, it is not 
yet possible from the observations to defini- 
tively exclude the possibility that WR 1 , WR 6 
and WR134 are binaries, and that WR137 has 
a third low-mass companion. However, the 
only viable binary scenario involves the for- 
mation of massive binaries with a high initial- 
mass ratio (up to 10), which seems to be ex- 
tremely unlike ly both observationally and the- 
oretically (e.g. Garmanv. Conti fc Massev _ 1980; 



iKobulnickv. Fryer fc Kiminki ll2007HBonnell 112007)) . 



5.2. The Non-Spherically Symmetric Wind 
Scenario 

If WR 1, WR 6, WR 134 and the WR star in the 
WR137 binary system are single stars, the peri- 
odic variability very likely originates in an asym- 
metry in their wind modulated by the stellar ro- 
tation. Spectropolarimetr ic observations of WR 6 
WR134 and WR137 (ISchulte-Ladbeck et al 



19911 : iRobert et al. I Il992t lHarries et al. I 11998) 



tering, indicating that the winds of these stars 
are not spherically symmetric. This is revealed 
by the depolarization of emission lines compared 
to the neighboring continuum. It was also found 
that the degree of depolarisation of spectral lines 
increases with decreasing ionisation state. This is 
interpreted to be a consequence of the ionisation 
stratification of hot stellar winds; as scattering 
lines with a higher degree of ionization are formed 
deeper in the wind, their polarisation level is closer 
to that of the continuum because of the higher 
density close to the core and therefore their level 
of depolarisation is lower. If a line has a strong re- 
combination component, wich is not polarised, its 
polarization level will be even lower and the depo- 
larization with repect to the continuum stronger. 
Finally, the level of linear polarization of the con- 
tinuum as measured in braodband light of WR6, 
WR 13 4 and WR 137 is extremely variable and pe - 



riodic (|Drissen et al. Ill989l; iMoffat et al. I [19931 



have shown the presence of an intrinsic contin- 
uum polarization component due to electron scat- 



This suggests that the wind density has a varying 
asymmetric distribution, such as density struc- 
tures that extend rather far in the wind. Unfortu- 
nately, there is no published high signal-to-noise 
spectropolarimetric nor continuum polarization 
observation of WR1. The only obse rvation so 
far was carried out by ISchmidt" ( 1988 ) and does 
not show a depolarization in its emission lines 
compared to the underlying continuum. The lat- 
ter, however, shows a significantly high level of 
polarization, although it remains to be demon- 
strated that this is intrinsic to the star rather 
than interstellar in origin. If the polarization is 
proved to be interstellar, this would explain why 
no depolarization is observed in the lines since 
the continuum is not polarized to start with. If, 
on the other hand, the light from the star is con- 
firmed to be highly polarized, a possible explana- 
tion for the lack of depolarization in the lines is 
that for some yet unknown reason, the region in 
which the lines arise shows a degree of polariza- 
tion as high as that in which the continuum arises. 
Also, the epoch- dependent nature of the changes 
could mean that relatively quiet periods can exist. 
Therefore, more spectropolarimetric observations 
are needed to verify wether the polarization of the 
emission lines of WR 1 varies with time. 

The possible existence of large-scale den- 
sity struc tures, s u ch as CIRs, was first pro- 



posed by iMullan I (|1984[ ). Following this idea, 
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Cranmer fc Owocki (1996) modeled the propa- 
gation of CIRs in a hot, radiatively, line-driven 
stellar wind. In that context, CIRs are caused by 
perturbations at the base of the wind, which in 
turn could be caused for example by a magnetic 
field or pulsations. These perturbations propagate 
through the wind while being carried around by 
rotation. This generates spiral-like structures in 
the density distribution that can lead to a char- 
acteristic, large-scale, periodic variability pattern 
in WR-wind emission lines, extremely similar to 
what is observed in the WR stars we discuss here 
( Dessart fc Chesneaull2002[ ). Thus, taking into ac- 
count the photometric and spectroscopic periodic 
variability, the spectropolarimetric observations 
and the low soft and hard X-ray fluxes, we con- 
clude that CIRs constitute an extremely likely 
interpretation. As for the epoch-dependency, it 
can easily be explained by the variable behaviour 
of perturbations that generate CIRs together with 
a finite lifetime of the structures in the wind. In- 
deed, even if the period of the variability caused 
by the motion of CIRs is always the same, the 
number and position of the CIRs may change de- 
pending on conditions at the surface. The origine 
of CIRs is still debated and no clear theoritical 
predictions of the typical lifetime of a CIR have 
been made yet. However, one can speculate that 
it will depend on the lifetime of the perturbation 
at the base of the wind and, to a lesser extent, 
on the flow speed in the wind. Observationaly, it 
could be determined from continuous observations 
during a great number of rotational periods (e.g. 
several weeks for WR6 and WR 134). 

5.3. Remarks on the putative CIRs in the 
wind of WR 1 

We believe that the best scenario to explain the 
photometric and spectroscopic variability we have 
detected in WR 1 is the presence in its wind of 
large-scale structures, most likely CIRs. In that 
case, each CIR would translate into a bump in the 
light-curve and a bump over the spectral emission 
lines. 

In theory, CIRs do not suffer from differential 
rotation and if the recurrence of the changes is 
caused solely by their rotation, they provide a di- 
rect measurement of the rotation period of the un- 
derlying star at the position in which they origi- 
nate, most likely close to the stellar surface. Thus, 



assuming a value for the radius of WR 1 of 2.2 R 
(a radius that c orresponds to a Rosseland opti- 
cal d epth of 20; Hamann. Grafener fc Liermann I 



2006), we obtain an equatorial rotational velocity 
of 6.5 km s . This is an order of magnitude lower 
than the values obtained for WR6 (40 km s _1 ) 
and WR 134 (70 km s" 1 ) when assuming a radius 
of R „, ~ 3i?0 (jHamann. Grafener fc Liermann 



2006). However, all these rotational velocities are 



in agreement with the very small values predicted 
for WR stars by m assive-star evolutionary m odels 
at solar metalicity (jMevnet fc Maederll2003l ). The 
order of magnitude difference between the rotation 
velocity of WR 1 and that of WR 6 and WR 134 is 
not completely surprising since the observed rota- 
tional velocity of a massive star at the WR evolu- 
tionary stage depends on several parameters such 
as the initial stellar mass, the mass-loss rate and 
the age of the star (time spent as a WR star). 
Interestingly, adopting a radius of 4.5 R© (cor- 
responding to the radi us of the hydrostatic core; 
iNugis fc Lamers" l l2000h for WR 137, we deduce an 



equatorial rotational velocity of 275 km s , i.e. 
more than half the breakup velocity. If confirmed, 
such a fast rotation for a WC star would render 
WR137 an interesting candidate for an eventual 
long-term gamma-ray burst. 

Ins pired by the analysis of Dessart fc Chesneaul 
(|2002l ). we can estimate the inclination of the 
CIRs in the wind of WR 1. Indeed, the maximum 
Doppler velocity that a bump associated with a 
CIR reaches during its motion on an emission line 
is v m ax = ivif r cos(9), where vif r is the velocity 
of the wind at the radii corresponding to the for- 
mation region of the observed emission line over 
which the bump is observed, and 6 is the incli- 
nation angle of the CIR with the line-of-sight. In 
Figure [T3l we track two bumps that reach v max i ~ 
-1300 km s _1 and v maX 2 ~ +700 km s _1 on the 
HeiiA5411 line. Assuming that the velocity law of 
the wind can be described by a /3-law for which 
the velocit y as a function of radius can be written 
as follows ( Castor fc Lamers~lll979l ): 



v(r) 



i?* 

r 



(4) 



where is the terminal velocity of the wind, 
the value of vif r for HeiiA5411 can be esti- 
mated using the emissi vity function described by 
Lepine fc Moffat I (|l999l) . However, that function 
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depends on the value of (3. Unfortunately, (3 can- 
not be estimated observationally for WR1, since 
no clumps are observed in the spectra. We can, 
however, give estimates of as a function of (3. If 
(3=1 , 2 or 3 the first CIR has an inclination angle 
of 0i =±40°, ±38° or ±30°, respectively, and the 
second bump has 6*2=±65°, ±64° or ±62°, respec- 
tively. When (3 is greater than 3, the velocity of 
the HeiiA5411 line formation region is lower than 
1300 km s -1 . Hence, under our first assumptions, 
we deduce that the (3 value for the wind of WR 1 
should be lower than 3. 

Of course, the association of the period of the 
CIRs with the rotation period of the underlying 
star remains to be confirmed. In the wind of OB 
stars, the presence of CIRs is thought to be at 
the origin of periodically recuring Discrete Ab- 
sorption Components (DA Cs) in P Cygni profile o f 



UV resonance lines (e.g. iFullerton et al. 1997 ) 



In many cases, the period of the outward moving 
DACs is found to correspond to v eq sini, where 
v eq is the rotational velocity at the equator and 
i the inclin aison of the rotational axis with the 
line of sight (Irlenrichs. Kaper fe Zwarthoed Ill988l 



Prini a 1119881) . But in the case of HD 64760, 



B0. 5 lb star, the CIRs rotate more s lowly than 
the stellar surface. Lobel fe Blommd (|2008h car- 
ried out a hydrodyna mical simulation of the CIR s 
for that star, following lCranmer fe Owockil(|l996f ). 
with the difference that they allowed the "spots" 
at the origine of the CIRs to move on the stellar 
surface. They were able to obtain a fairly good fit 
to the lpv of the SiivAA1394,1403 and concluded 
that, for that star, the origin of the CIRs must be 
the interference pattern of a number of non-radial 
pulsations at the surface of the star. 

So far, there is no model for the "spots" that 
are at the origin of the creation of the CIRs and 
there are only a limited number of datasets avail- 
able that could potentially be used to confirm 
either the pulsational or the magnetic origin of 
CIRs. The range of frequency expected for pul- 
sations in WR stars is not well determined. Up 
to date, only one detection of a relatively sta- 
ble 9.8 h pe riod that can be attribu t ed ei ther 



to g-modes (ITownsend fe M acDonald 200 6J) or 



stange-mode pulsations (Dorfi, Gautschv fc Saio 



2006) has been claimed in WR 123 ([Lefevre et al. 



a magnetic field has been achieved in a WR star 
so far; only an upp er limit of ~ 25 Gaus s for WR 6 
has been claimed ( St-Louis et al.l[2008h . 



2005bh . No pulsational period of a few days is cur- 
rently known nor predicted. Also, no detection of 
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